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SUMMARY 


rrA 


The Lindstedt perturbation method is applied to the motion of an artificial Earth satellite that moves 
along a geodesic of the Schwarzschild metric of general relativity theory. It is shown that both the ampli- 
tude and the frequency of the first-order approximate solution obtained are affected by the nonlinearity of the 
relativistic term appearing in the equation of motion. The approximate periodic solution is compared with 
the solution for the motion of a synchronous artificial Earth satellite (ATS 3) in a Newtonian force field. 
Under the assumption that the coordinate time and the initial conditions are the same in both systems, it is 
deduced that the maximum of the magnitude difference between the two calculated radii is of order 1.6 cm 
after one-half orbit and is, therefore, too small to be detected. 
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INTEGRATION OF THE REUTIVISTIC EQUATIONS 
OF MOTION OF AN ARTIFICIAL EARTH SATELLITE 




by 

Abolghassem Ghaffari 
Goddard Space Flight Center 


INTRODUCTION 

In the theory of general relativity, the external gravitational field of a spherically symmetric massive 
body M whose center lies at r = 0 is represented by the static Schwarzschild metric 

ds^ - ~ ~ ^ r^{dO^ + sin^^d<?S^)j . (1) 

In this metric, an exact solution for a set of discrete bodies is possible only when one of the bodies is of 
finite mass and all of the rest are of infinitesimally small mass (Reference 1). If the massive body is taken 
to be the Earth, then an artificial Earth satellite may be considered as an infinitesimal test particle whose 
gravitational field may be neglected. 

The complete solution for the relativistic effects of the combined fields of the Earth and the Sun acting 
on an artificial Earth satellite is complex. The metric must involve both the principal masses, and since 
the Earth is orbiting the Sun, the field is no longer static. It is well known (Reference 2) that, in general, 
the motion of an Earth satellite can be derived with a great deal of precision through consideration of the 
effects of the Earth and the Sun separately. The relativistic effects of solar gravitation have been calcu- 
lated by Papapetrou (Reference 3) and Corinaldesi and Papapetrou (Reference 4). The use of artificial 
Earth satellites seems more suitable than the use of solar satellites for testing the effects of general rela- 
tivity because the Sun has many planets and their mutual attractions are significant. 

For convenience, it is assumed that the Earth's gravitational field is spherically symmetric and that 
the Sun's gravitational effect upon an artificial Earth satellite is negligible. With these assumptions, one 
knows that the motion of an artificial Earth satellite is given by the equations of the ordinary geodesics of 
the Schwarzschild metric (References 1 and 5) which combined into one simplified form is 


*Some of the results of this paper have been presented at the Sixth Semiannual Astrodynamic Conference held at 
Goddard Space Flight Center, Greenbelt, Maryland, on November 7 and 8, 1967. 
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or 



( 2 ) 


where 



and 

r = distance of the satellite from the center of the Earth, 

Af = mass of the Earth, 

V = speed of the satellite, 

c = local speed of light, 

/i, /3 = constants of integration, 

<56 = true anomaly, 
and 

G = gravitational constant. 


INTEGRATION OF THE EQUATIONS OF MOTION 

The classical equation [equation (2)] has been investigated and discussed completely in most standard 
texts on celestial mechanics and general relativity (References 5, 6, 7, 8, and 9). The rigorous integration 
of equation (2) leads to the Weirstrass elliptic function a (^, gg, gg), which satisfies the differential 
equation 

where gg and gg are constants. However, in practice, an approximate integration of equation (2) gives the 
advance of the perihelion. The solution of equation (2), with the Sun as the central body and a planet as 
the test particle, gave the perihelion advances in 100 years for the planets Mercury (43'.'03), Venus (8'.'64), 
Earth (3'.'84), and Mars (1'.'35) predicted by Einstein's law of gravitation. A complete discussion of equation 
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(2) has been given by Chaz;y (Reference 6), and exhaustive research of the orbits defined by equation (2) 
with the location of their singularities has been performed by Hagihara (Reference 10). McVittie (Reference 
5) discussed and solved equation (2) by a different method. 

Differentiating equation (2) with respect to the true anomaly <f> and setting mc^/h^ = 1/p, gives 

H- u =— + Smu ^ , (3) 

P 

provided du/d^ has only isolated zeros (i.e,, eliminating the circular-orbit solutions). In the Newtonian law 
of gravitation, the right-hand sides of equations (2) and (3) are quadratic and linear in u, respectively. The 
presence of corrective terms 2mu^ in equation (2) and in equation (3) arises from Einstein's law of 
gravitation. 

Equation (3) was first derived by Eddington (Reference 7). Using a method of successive approximation, 
he obtained a Keplerian orbit as a first approximation by ignoring the small term 3mu^. Then, substituting 
the first approximation into the small term he arrived at a second approximation with a secular term 

describing the resonant case. Bergmann (Reference 9) also considered equation (3) and, applying a Fourier 
series procedure, he obtained the perihelion advances of the planets without giving the explicit expressions 
for the planets’ orbits. 

In the sections that follow, 

(1) Equation (3) is integrated by the classical Lindstedt method to give a first approximate periodic 
solution starting from the perigee. 

(2) The approximate periodic solution is compared with the solution obtained for the motion of an arti- 
ficial Earth satellite in a Newtonian force field. 

(3) A numerical estimate of the comparison is made with the assumption that the central body is the 
Earth and the test particle is a long-life, synchronous Earth satellite such as the Application Technology 
Satellite 3 (ATS 3). 

The advantages of equation (3) over equation (2) are twofold: 

(1) Though equation (3) is of second order, it is, nevertheless, nonlinear in u only, whereas equation 
(2) is nonlinear in both u and du/d(f>, and its integration leads to the Weierstrass elliptic functions. 

(2) \Ve are primarily concerned with the periodic solution and its behavior of periodicity. The left-hand 
side of equation (3) clearly suggests the application of classical and modern approximation methods for the 
derivation of approximations to periodic solutions. 

APPLICATION OF THE LINDSTEDT METHOD 

Let us assume that we wish to find an approximate periodic solution of equation (3) satisfying the fol- 
lowing initial conditions (i.e., starting from the perigee): 
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In order to find such a solution and to investigate the perturbation in the basic frequency arising from the 
presence of the relativistic term on the right-hand side of equation (3), and also to eliminate the secular 
term, the Lindstedt method (Reference 11) is applied to equation (3). 

The Lindstedt method consists primarily of a change of the independent variable from ^ to another in- 
dependent variable a such that the determination of the available unknown coefficients enables us to elimi- 
nate gradually the secular terms in the subsequent approximations (Reference 12). As the frequency is 
altered, it is of advantage to replace the independent variable ^ with a new independent variable a in a 
manner defined by the relation 

- a(l + c^€ + Cge^ (5) 

where c^, Cg, . . . are unknown coefficients and is an arbitrarily small positive parameter. The smallness 
of the gravitational radius m enables us to suppose that e = 3m. Then, equation (3) becomes 

^+(l + Cj,+C^e^ + ...f (u-^-eU^y-0, 
or 

■^-^+u -—+((- — ^ + 2c- u - u^) + o(f^) = 0 . (6) 

da^P\P^/ 

The solution of equation (6) can be written as a power series with respect to the small parameter e : 


u(e, a) 


00 

n=0 


( 7 ) 


To limit ourselves to the first-order approximation, the series is written 

u = Uq + e-Uj + (8) 

Using equations (6) and (7), one finds that the leading term Uq is a solution of the unperturbed equation 

d^Ur 


da^ 


0 . 1 
+ Un = — 


( 9 ) 


and that satisfies 


d^u^ 

da^ 




2 + 
2p^ 


2e 

— + - 5(1 - c.p) cos a + — ^ cos 2a. 

P p^ ^ 2p^ 


(10) 


The unknown coefficient c is to be determined such that no secular term appears in the solution of equation 
(10). Hence, we choose so that 


c 


1 “ 




( 11 ) 
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and, therefore, equation (10) assumes the form 


d^u 

da^ 


u, = - 2 + cos 2a) . 

^ 2p‘^ 


( 12 ) 


Equation (12) has the solution 




satisfying the initial conditions: 


«l = 


e^ — 2 

5^(1 - COS 2a) 

2p‘^ 


(13) 



(14) 


Therefore, the first-order approximate periodic solution of equation (3) satisfying the initial conditions 
[equations (4)] is 


where 


, , 2\ 1 + e cos a , 3m(e^ “ 2) , , g, 

u = Uf, + fu, + o{e'^) = + i — 5 — i(l - cos 2a) + o(e‘^) , 

“ P 2p‘^ 

= o(.")=(l-»£V+ 

1 + Cjf + . . . 1 + V p 


(15) 

(16) 


We can deduce that the nonlinearity shown by the relativistic term in Equation 3 affects not only the ampli- 
tude but also the frequency of the solution. 

Equations (15) and (16) show that the change in the frequency depends upon the amplitude a and the ec- 
centricity e of the Keplerian orbit and also upon the parameter e = 3m , a property of the periodic solutions of 
all nonlinear autonomous differential equations. The period of the approximate periodic solution of equation 
(9) is 2 t 7, i.e., the orbits are closed. The period of the exact solution of equation (3) differs from by a 
small angle d which is the difference between the angular positions of two succeeding perigees and is given 
by 

5^ , 2n = 2n(l +^)- 2.=^= (17) 

1 - (3m/p) \ P J P a(l - e2)\rev/ 

Therefore, the precession of the perigee of a satellite orbit obtained by this method amounts to Qinn/p 
rad/rev, which is in close agreement with the precession predicted by Einstein’s theory of general relativity 
as well as with the gravitational theories of Whitehead (Reference 13) and Birkhoff (Reference 14). 
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COMPARISON WITH NEWTONIAN FORCE FIELD 
AND NUMERICAL ANALYSIS 


Now let a second artificial Earth satellite with the same characteristics as the first move in a planar 
elliptic orbit according to Newtonian laws of motion. If and 0 are the classical polar coordinates of the 
second satellite in the orbital plane and = 1/r^, the classical Newtonian equation of motion is 




1 + e cos cj^ 


(18) 


where the true anomaly 9 S is measured from the perigee and p = a(l - e^) is the semi-latus rectum; a and e 
are the semi-major axis and the eccentricity of the Newtonian orbit, respectively. We assume that both sat- 
ellites start from the perigee, i.e., they both satisfy the initial conditions: 






('"w)o“Cr)o' 


(“w)o - (^ r)o ’ 


(19) 


(“n)o = (“r) 


0 ’ 


where and stand for relativistic and Newtonian solutions, respectively. 

For a comparison of the approximate periodic relativistic solution [equation (15)] with the classical 
Newtonian solution [equation (18)] , the difference Au is formed: 


Au 




e r / /. 3/n\ ,~] 3m(e^ - 2) 3m\ ^ 

: — cos 0 ~ cos ( 1 )o ^ — n 1 - cos 2(1 )(,6 

PI \P/J 2p^ I \P/J 


( 20 ) 


The difference Au is a function of the Keplerian orbital elements and the gravitational radius m of the cen- 
tral body. Since the central body is taken to be the Earth, m ~ 0.443 cm, and the parameters a and e are 
constant for one revolution. Therefore, for a single revolution, the difference Au = Au(<; 6 ) is a function of 
the true anomaly 0 only. 

For a numerical estimate of Au(^), equation (20) was applied to orbital data of the ATS 3, which had the 
the elements 


and 


a = 6.6109161 x R , 
e = 0.0001703684, 


where R = 6.371 x 10® cm is the mean radius of the Earth. Setting 


~ ^ ~ 




( 21 ) 


where a denotes the difference between the relativistic and the Newtonian radii, one finds that 
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( 22 ) 


Au ~ u 


N 



a 



The results are listed in Table 1 and plotted in Figure 1 for the range 0 < 0 < 190° The function 
being a continuous function of begins to increase again after 180°. Table 1 and Figure 1 show 
that, for the case of the ATS 3, the relativistic correction to Equation 3 is too small to be detected, and the 
maximum value of Au ~ 1.5 x 10“^^ is obtained for = 89°59‘59'.'148792708 (correct to the nine decimal 
places given in seconds). Therefore, the maximum radial deviation is of the order 

r^(1.5 X 10“^®) = (4 x 10^ cra)^(1.5 x 10~^^) 1.6 cm, (23) 


which is too small to be detected. However, it is important to separate the relativistic effects from those 
attributable to other causes, such as the oblateness of the Earth, magnetic fields, the influence of high- 
altitude winds, and so forth. 


Table 1-Variation of Au(0) in terms of the true anomaly 


(deg) 

Au(0) 

0 

0 

10 

.4517717214 x lO"^^ 

20 

.1752594324 x 10“^® 

30 

.3745562899 x 10~^® 

40 

.6190285723 x 10”^^ 

50 

.8791879840 x 10~^^ 

60 

.1123653936 x 10”^® 

70 

.1322938634 x 10'^® 

80 

.1453003889 x 10”^® 

89 59'59'.'148792708 

.149816062806 x 10"^^ 

90 

.1498160626 x 10"^® 

100 

.1452961373 X 10" 

110 

.1322857458 x 10" 

120 

.1123541627 X 10" 

130 

.8790553806 x 10"^^ 

140 

.6188892964 x 10" 

150 

.3744260732 x 10"^^ 

160 

.1751553208 X 10“ 

170 

.4511663842 x 10"^° 

180 

-.1050408293 x 10”27 

190 

.4525317707 x 10“20 
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In a later paper we expect to report on the application of the averaging method of Bogoliubov and 
Mitropolsky (Reference 15) to equation (3), which leads to a different approximate solution, and its compar- 
ison with the solution obtained here [equation (15)]. A comparison of the obtained solution with that of 
Lass and Solloway (Reference 16) will also be included. 
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